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Abstract 

A new mechanism different from the spin accumulation picture is proposed for the current 
induced magnetization switching in magnetic tunnel junctions by taking into account the effect 
of the electron electron interaction. We found in tunnel structures the possibility of an enhanced 
spin switching effect that, when normalized with respect to the current, is much bigger than that 
in multilayers. Some recent experimental results show evidence for the present picture. 
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Spin polarized transport such as the giant magnetoresistance (GMR) and the tunnelling 
magnetoresistance (TMR) has recently received much interest Q |], (7|. This 

is partly motivated by technological applications such as magnetic sensors, hard disk read 
heads, magnetic non-volatile random access memories (MRAM) and the possibility of inte- 
grating magnetic elements into semiconductor microelectronics circuits. The crucial physics 
that goes into the design of these structures is that of the transport of electrons over dis- 
tances short enough that the spin memory has not been lost. One of the effects that is 
manifested is that of spin induced switching of magnetization. 

When a polarized electric current goes from one magnetic element to the next, the mag- 
netization carried by it when it leaves the first element will exert a torque and rotate the 
magnetization of the second element. This is the spin accumulation picture for magneti- 
zation switching. Following earlier theoretical suggestions |[ | a number of experiments 
[|. |5| have recently demonstrated this effect under high current densities. If the current 
required can be reduced this rotation may provide a way to write a bit in a memory device. 
In this Letter, a different mechanism is proposed for magnetization switching in magnetic 
tunnel junction which is expected to be realizable under much lower current densities. 

Intuitively, for tunnel junctions, the resistance is higher, the current and the magnetiza- 
tion carried by it is reduced. According to the spin accumulation picture the magnetization 
switching effect will be reduced. However, there are additional physics that needs to be 
included, which leads to a new mechanism of the magnetization switching effect. The tun- 
nel junction can be thought of as a capacitor with charges localized at the metal-insulator 
interfaces within a screening length A. To consider the capacitance correctly it is essential 
to include the effect of the electron-electron (e-e) interaction in the nonequilibrium transport 
of the system. We found that very interesting physics happens after taking into account the 
e-e interaction. There are two very different length scales in the problem: The bare spin 
diffusion length l s f, of the order of 100 A or more, is much larger than the bare screening 
length Ao which is of the order of an A. Under steady state non-equilibrium conditions, 
the external voltage induced a non-uniformed charge dipole layer at the metal-insulator in- 
terfaces. These charges are sums of two terms. In addition to one decaying with a length 
scale of the order of the screening length that one encounters in the equilibrium situation, 
there is another one decaying with a length scale of the order of the spin diffusion length. 
The magnitude of the density of the latter is smaller than that of the first but they con- 
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tribute equally to the screening electric field due to the long range nature of the Coulomb 
potential. In addition, there is a magnetization dipole layer also induced at the 
interfaces. The magnitude of this magnetization dipole layer is much larger than that of 
the charge dipole layer by a factor of l s f/Xo- This magnetization dipole layer is zero under 
equilibrium conditions and is nonzero only under steady state conditions. It generates an 
effective switching field. In particular, the magnitude of the induced magnetization is con- 
trolled by the voltage, not the current, implying a possibility of achieving an enhanced 
switching effect in tunnel structures that, when normalized with respect to the current, 
is much bigger than that caused by the spin accumulation picture. We thus expect this to 
have a strong effect in tunnel junctions where the voltage involved is higher while the current 
density remains low. There is some experimental evidence for the present mechanism, 
although definitive studies remain to be carried out. We now describe our results in detail. 

The system we have in mind is a ferromagnetic tunnel junction where the two interfaces 
between the ferromagnet-insulator-ferromagnet sandwich structure are assumed to be at 
z = id/2. We assume the z axis to be perpendicular to the faces of the tunnel junction. 
The initial magnetization are assumed to be in the x-y plane with an orientation given by 
= cos $ e x + sin $ e y for the ferromagnet on the right hand side and pg = e x for the 
ferromagnet on the left hand side of the sandwitch structure. Our strategy is to assume 
magnetization uniform in the x-y plane on both the left and the right, consider a current 
going through this junction and calculate the magnetization change <5M under steady state 
nonequilibrium conditions. The effective switching field is then estimated by as H s = 5~M./x 
where \ is the magnetic susceptibility. Because the work functions of the ferromagnets on 
opposite sides of the junction may not be equal, at zero external bias there will be a charge 
dipole layer formed at the interfaces. What we are calculating here are the changes from the 
zero bias situation. The experimental structures usually possess edge domains where the 
switching starts. The magnetization is thus not completely uniform in the x-y plane. To 
bring out the essential physics, we shall not consider this complication in the present paper 
but we hope to come back to this in the future. 

Our starting point is the equation of motion of the charge and the magnetization. For 
the charge, it is just the equation of charge current conservation 
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where J e is the total current. For the magnetization M, the equation takes the form of the 
classical Landau-Lifshitz equation 

— - 7 MxH + V-J M = — (2) 

where 7 is the gyromagnetic ratio, Jm is the spin current (tensor), and H is the effective field 
describing the precession of the magnetic moments given by H = H c + H an + H dip + H ex . 
H cx = JV 2 M is the effective field due to direct exchange; H an = fTM denotes the influence 
of anisotropy energy; H e represents the external field; and Hdi P denotes the dipole-dipole 
interaction, r is the relaxation time, describing the relaxation of the system towards its 
local equilibrium value of magnetization. Eq. (1) and (2) involves magnetization and charge 
currents, to which we turn our attention. 

To derive the transport equation for the vector magnetization, we consider a local coor- 
dinate system so that the electrons are quantized in the y' direction. We write down the 
transport equations in this local coordinate system and then rephrase it in covariant form, 
thus obtaining the general transport equation. According to Fick's law, the current for spin 
s is given by 

-eVV + sfiB^Hy, -D s Vn s , (3) 



_ n s /j, s 

<> s — 

e 

where \ib denotes the Bohr magneton, n s , D s , and /i s are the number density of charge 
carriers, the diffusion coefficient, and the electron mobility, respectively, e < is the electron 
charge, s — ± denote different spin orientations, and V = V e + W, with V e the electric 
potential describing the external electric field and W the local electric (screening) potential 
due to the other electric charges determined self-consistent ly by 

W(r) = J d 3 r'U{r - r')dn(r') (4) 

with U the Coulomb potential. The total number density of charge carriers and y' component 
of magnetization are given by n = n + + n_, M y i = /is(n + — n_). 

From (3), we obtain expressions of the total charge current J e and magnetization current 
Jm- They are, after linearization and rephrasing in covariant form, given by 

J e = -aVV - a ie £DV5n - — /3fDV(<M • p ) 

' (5) 

Jm = -^mV(Vpo) - ££>V(M - a 2 (3^DV{6np ) 
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where ff = eV n s fi s , om = /4b 52 sn s^s, P = — — , Po — tttt with M the local equilibrium 
s s Vbo- |M | 

2 

magnetization, and £ = with D = | D s and y the Pauli susceptibility at the Fermi 

e 2 xD 2 V 

energy, cti and «2 are two renormalized phenomenological parameters, with a\ ~ «2 ~ 1 in 
the limit of strong ferromagnet. 

Substituting the expression for Jm into the Landau-Lifshitz equation (2) we obtain the 
relaxation equation for M: 



V 2 cM - W SNL + CPo x (V 2 <M - W SM) = -(3fi B p (a 3 V 2 Sn - (6) 

's/ ^0 



where only H cx = JV 2 M = JV 2 5M and H an = ifM are kept in the precession term 

e 
— ( 
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7M x H, and use has been made of Gauss' law: V 2 V = V 2 W = Sn. The bare spin 



diffusion length l s f and the bare screening length An are given by 

Z 2 = r^D and A 2 = ^ 
J a 

respectively. Other dimension^ parameter are C = '^f and « = M while a 3 is 
again a renormalized phenomenological parameter collecting all V 2 5n dependence including 
the change of D (such as is caused by the dipole layer at zero bias due to the difference of 
the work functions of the material on the left and the right) near the interface. 
The charge current conservation (1) yields, for the steady state, 

Sn - ai fi B V 2 Sn - (3V 2 (6M ■ p ) = 0, (7) 

which, together with Eq.(|), describes the distribution of the charge and magnetization 
on opposites sides away from the tunnel junction in terms of their values at the junction. 
The values of the charge and magnetization densities at the junction can be determined by 
matching boundary conditions across the barrier. We first solve these equations in the metal 
part of the junction. These solutions determine the charge and magnetization dipole layers. 

We expect the charge and magnetization dipole layers to decay away from the interface 
with length scales controlled by the spin diffusion length and the screening length. Because 
of the vector nature of the magnetization, there are three normal modes by which they 
can decay away from the interface. Including the charge degree of freedom, there are four 
normal modes that one can consider. For ferromagnetic metal on the right hand side, we 
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thus consider the following ansatz: 



5n R = J2 Sn?oe- {z ^ )/k , 5M R = J2 <Mf e~M)A 



(8) 



i=i 



where the superscript R denotes the right hand side. Inserting these into and letting 
the coefficients before the exponential scaling functions vanish for steady-state solutions, we 
get the renormalized decay lengths 



h 

with 



An, h 



P 2 hf, h = («+ + ia-)l 8 f, h = ll = (a+ - iaJ)l sf (9) 



a± 



:i+c 2 « a )(i+c 2 )±(i+«c 2 



2(l + C 2 /« 2 ) 

As we indicated in the introduction, the screening length and the spin diffusion length are 
renormalized. As we shall see below, £3 and £4 correspond the length scales with which 
the "precession" dies away from the interface. The charge densities can be related to the 
magnetization densities by 

A 2 , 



1-/3 2 5M R 



5n R 



6n R = 0. 



(10) 



P(ai - a 3 ) hb ' " 3) l\ (J>b 

Note that Sn R is much smaller than SM^ because Ao << h- Inserting the "eigen-solutions" 
into equations (||), we finally obtain analytic expressions for the dipole layers: 



5n R = 5nf e-^y h + Sn^e^'-^ 
(M R = v R 5M R e^ z -iV h +v R 5M R e-( z -^ 



(11) 
(12) 



where 



a_z 



77 = — , c = a, +a_, q = sin (p' sin $ e x — sin (p' cos $ e v + cos (p 1 e z , <p' = tp — . (13) 
cr c z l s f 

5Mtfi, with i = 1,2, 3, and <p R are to be determined later. Terms of the order (An/Z s /) 2 have 
been neglected since / 2 j >> Aq. As advertised, the charge dipole layer is the sum of two 
terms, one decaying with a length scale of the screening length; the other, the spin diffusion 
length. The vector magnetization dipole is now a sum of three terms. The first two ( $VI{^, 
SM. R q ) are along the direction of the original magnetization; the last one is perpendicular to 
the direction of the original magnetization and describes the precession of the magnetization 



around the original axis. Again, the first two terms correspond to decay lengths of the order 
of the spin diffusion length and the screening length, while the precession term only decays 
with a length scale of the order of the spin diffusion length. 

With Sn L , SM. L , Jg , and obtained similarly for the ferromagnet on the left hand side, 
we can next focus on the boundary condition across the insulator. The charge neutrality 
condition 



6n R dz + / 2 5n L dz = 



yields 



hi 



0. 



We have, for simplicity, assumed the identical parameters D R = D L 
a R = a L = a, etc. The continuity of charge current J e and the spin current Jm implies 



(14) 
(15) 



z=-d/2 



3 R 



z=d/2 



M 



z=-d/2 



1 R 



z=d/2 



(16) 



The charge and magnetization currents (||) can be expressed in terms of the SM^, etc as: 



iR 



trE, 



ext 



3 R = aE ext ^/3 P * + K D "f ~ as) e zP R 5n R e-^ 



+ (1 f3 2 )ZD ezp R SM R e -(z-l )/h + iR GzCl R 5M R e -v( Z -i)/is f 
h ' cl s f 

where E ext = E ext e z is the external electric field inside the conductor, and 



(17) 



(18) 



sin cp" sin $ e x — sin cp" cos $ + cos cp" e z , cp" = cp' — sin 1 



(19) 



Note that q m ■ po = ■ Po = 0. The currents J e and Jm are assumed to satisfy the boundary 
conditions 
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(G R 5n R - G L s Sn 



erd 



(G R SM R ■ p R - G L D m L ■ p£ 



erd/j B 



^■J°m = -^-iG R 5a R V R - G L D 5n L p%) - —^—(G R SM. R - G L S <M L ) 



(20) 
(21) 



e 2 rd K u ^ u ^' e 2 rd' 

where G R and G R (Gj) and Gg) are phenomenological parameters describing combinations 
of inverse densities of states of the ferromagnet on the right (left) hand side ||, and r is the 
resistivity of the junction (insulator). The AW terms are much smaller and neglected for 
the system. 



^From equations ([15]), (00), and (gO^gl]) we can solve for 5nf , 5n? , 5M^ 0: 5M 3 L , 



Since 



7 \ r^ L ( R ) 



<< l "v^ << 1, 

Z2 l s f ' -De 2 r rcr 

terms of higher order are neglected. We find that the dominating term 5M 2 o is proportional 
to the voltage rdJ e . The resistance of the junction, and not that of the metal enters into 
consideration. More specifically 

5M£ = -6M* = -^4- (22) 

All other terms such as SM^ , SM^ Q , and Sn^o are much smaller than SM^q- They are given 
by 

<M& = -5M« = = §6M» = Cl f n *r* 5M^ <f = <p* = cos- 

(23) 

where it has been assumed that Pg = e^, = cos $ e x + sin $ e r The charge densities are 
related to the magnetization densities through Eq.fllQl). 

Eq.(^2|) is our main result. It shows that an effective switching field H s = 5M^/x can 
be generated which is basically controlled by the voltage and, when normalized with respect 
to current, is much bigger than that caused by spin accumulation picture. By controlling 
the sign of J e , the sign of this field can be changed. In this calculation, we have assumed 
the magnetization is uniform in the x-y plane. In a small element, there are edge domains 
where the switching first starts. We expect that this effective field will also act on the edge 
domains. 

There is some experimental evidence of a much stronger current induced switching effect 
than expected based on the spin accumulation picture. Kohlstedt and coworkers have 
recently studied experimentally a three terminal multiple junction Fx — I\ — F2 — I2 — F3 
with leads at Fx, F 2 and F 3 . They found that the magnetoresistance for F 2 — I 2 ~ F 3 is 
changed as a current is passed through Fx — I\ — F%. They interpret the change in the 
magnetoresistance ratio in the second tunnel junction as a change in the magnetization in 
F2 when a current is passed through the first tunnel junction. They found that the change 
in magnetization in F2 is much bigger than that caused by the spin accumulation 
effect! Furthermore, the change in the resistance of the second junction is increased as 
the resistance of the first junction is increased, opposite to the trend according to the spin 
accumulation picture. Their result is consistent with the present picture. 
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In conclusion, We have studied the transport of the vector magnetization through a 
magnetic tunnel junction. Because of the induced charge and magnetization dipole layers at 
the interfaces, we found in tunnel structures the possibility of an enhanced spin switching 
effect that, when normalized with respect to the current, is much bigger than that caused 
by the spin accumulation picture. 
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